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INTRODUCTION 


It is universally assumed in the text-books on 
suspension bridges that the shape of the curve 
assumed by the cable isa parabola on the suppo- 
sition that the weight of the structure and the loads 
are uniformly distributed in a horizontal direction. 
Neglecting for a moment the moving loads, the 
weight of the structure is made up of the weight of 
the cable and the suspenders, and the weight of the 
roadway and the stiffening strusses, it being assumed 
that the stiffening is not of the braced-chain type. 
We could simply classify this as the weight of the 
cable and the weight of the roadway, the former 
being uniformly distributed along its length and the 
latter horizontally. It is well-known in statics that 
if the former be inappreciable and can be neglected 
the curve of equilibrium is a parabola; while if 
the latter be neglected in comparison with the 
cable weight the curve is a catenary. It is hardly 
necessary to point out that the curve will still 
be a parabola even if the former be not neglected 
but taken to act horizontally. I propose to examine 
in this paper how far this assumption is justifiable 








* Paper presented at the 26th Annual Session of the Mysore 
Engineers Association held in March 1934. 
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and what is the error made by the assumption in 
calculations relating to suspension bridges. In view 
of the fact that every book, either on applied 
mechanics or suspension bridges, makes _ this 
assumption without any attempt at justification, it 
is worthwhile to investigate how far it is correct. 
Before proceeding to this discussion, it might be 
observed that greater errors are likely to result from 
this assumption, the greater the weight of the cable. 
The following *table giving a comparative statement 











Net weight of Gross fixed 
Name. | Span in ft. chain per ft. load per ft. 
run. of span. 
Clifton 702.25 | 1560 Shee 
Freiburg 870 167 760 
Niagara 821.3 | 820 2032 








for three suspension bridges shows that the cable 
weight is quite appreciable when compared with the 
gross fixed load and it may consequently be of some 
practical use to undertake this investigation. 
Further, it is also evident that large values of the 
ratio of sag to span introduce greater errors in the 
assumption and in practice this ratio varies from 
0.2 to 0.05 such that it might be as high as 1/5. 


§ 1. THE GENERAL DIFFERENTIAL EQUATION 


To determine the form of the cable for any 
loading, consider vertical loads applied on a cable 
suspended between two points. The cable will 
assume a definite polygonal form determined by the 


*“ Table from Rankine : Civil Engineering, p. 582, 
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relations between the loads. The end components will 
be inclined and will have horizontal components H 
Simple considerations of static equilibrium show 
that H will be the same for both end reactions and 
will also equal the horizontal component of the 
tension in the cable at any point. H is called the 
horizontal tension of the cable. It is easily shown 
that the cable.curve is the, bending moment diagram 
for the applied loads. If the loads are continuously 
distributed the funicular polygon becomes a conti- 
nuous curve and if w be the load per horizontal 
linear unit-at any point having the abscissa x, we 
can obtain the equation for the vertical shear at any 
section x of the span and on differentiating this 
equation we obtain, 

dx 


=a 
and on taking the y-axis vertically upwards this 
becomes 
oe dy. _ ’ 

Hf dxtinae (1) 
which might be taken as the differential equation™ of 
the equilibrium curve in the general case. Let g. be 
the weight of the cable per unit ,length (supposed 
uniform) and w» the total uniform moving and dead 
loads per horizontal foot ; equation (1) then becomes 


2 
Ho, = g. secK + Wo (2) 


where « is the inclination of the cable to the horizon- 
tal at any point. Eqn. (2) is the general differential 
equation of the problem. If go be assumed to bg 








“ See for e.g., W. H. BurR: Suspension Bridges p. 61. 
Eqn. (13) or JOHNSON, BRYAN and TURNEAURE: Modern 
Framed Structures. Part II, pp. 187-88. We shall .use 
throughout the notation of the former book. 
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uniform horizontally instead of along the curve, we 
can reduce the equation to 


2 


= + 
H a Ww (w=wo ' So) 


and taking the initial conditions at the lowest point 
of the cable to be: 


x =0 dy =0, y = 0, this equation gives 
ne 
wx? _ wx 
Dy a (2a 


which is the standard equation to the parabolic cable. 
The tension T at any point is given by 


é‘ 
T = H Sec « = HA/, . (2) 


If the points of support are at the elevation y=f 
and the span = 7, H takes the value 
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we _ wh a 
H= Picts oe M1 +(¥) 
If L be the entire length of the cable we have 


= l/s (#2)? 4 by : ] * 
L. =2,) [ }+ Fe | axe (2u+sin h2u) 


ia font sit ddd lens 


Developing the radical] by the binomi 
three terms and performing the integr 


e Boda 34 £5) 

L if 1+8 rose fo | 

. Lhe equations for y, H, T, Lare the fundamental 
equations for the parabolic cable. 


al theorem to 
ation indicated 








“n = f/land u is defined by sin hu = 4n,’ 


This elegant 
form for L is to be found in Steinman: Suspen 


sion Bridges p. 5, 
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If in equation (2) we neglect w, we obtain 


- Acai ahaa 
HS = go sec & =.A/ ie (42)" 

This is the differential equation of a catenary and 

its solution gives 


evi Ce ox 
ae ao (i ay “) where c= go/H 


and we can derive corresponding equations for H, T 
and L as in the parabolic cable. 


I now consider the general differential equation 
(2) where neither w, nor go can be neglected (it 
must be observed that assuming a uniform horizon- 
tal distribution of g, is mathematically equivalent to 
neglecting it) and obtain the exact mathematical 
solution. In the next article, I derive an approximate 
form of the solution suitable for practical use and 
enabling us to compare the formule obtained for a 
strictly parabolic cable. 


We shall write (2) in the form 
d*y _ 


H qe sec © + wo (replacing * by © ) 
x 
r d*2= asec O+d. ( a=" ,b=*-) 
oO d x2 H H 


and assume, as always realised in practice, that b>a. 
Putting p=dy/dx, this further reduces to the first 
order differential equation: 


LI POLES 
42 — r/ (3) 


Since p = tan ©, dp = sec? ©. dO and (3) can be 
written as 


sec? @. dO _ Dep) Pe iA 25 
b + asec ree OF cos © (at+b cos O) - 





dO _ bd =—ad 
Oe teoeieiae a+bcos 8 a 


and integrating we have 





a* = log (sec © + tan O)- 
b+ o> b2—a? sinO 
a a i ae | a@ COS r/ a’ si (4) 
Vig — at ey 


the constant of integration being given by x = 0 for © = 0. 


a+bcos © 





Again oP a 1 + p* + bcan bewritten as 
x 


pdp tan © sec? Od O 
ay/i tpt sald 1.566 a+b 








= dy 


This can be easily integrated and the result is’ 


+ b 
a= —1) + a wrt | 
a? y = a(sec 0-1) + b log rin ma S 


taking y = 0 for 9 = 0. 


Equations (4) and (5) give the parametric equations 
of the curve of the cabie in terms of ©, but they are 
in quite an unmanageable form fot any practical 
purposes. Before deriving approximate equations 
from these, it is instructive to see how these reduce 
to the parabola and catenary for thecases a=O0 and 


b=0 respectively. 

Case b= 0: Equation (4) takes the simple form 
ax = log (sec 9 + tan @) 

and equation (5) reduces to 


ay = sec © —I1, since the argument of the 
logarithm in equation (5) can never become infinite. 
Eliminating © between these equations we have 
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= 1( ax —ax ) 
Y 2a\e +e —2 


which is the equation to the catenary 


Case a = 0: In this case equations (4) and (5) reduce to 
identities but by dividing them throughout by @ and 
a2 respectively and applying de L’Hospital’s Rule 
for the evaluation of indeterminate forms, they 
take the form 


tan © | 


— ole 


and » =>, tan? G 


bx? 
which give » =F , the equation to the parabola. 


§ 2. APPROXIMATE SOLUTION OF THE GENERAL 
DIFFERENTIAL EQUATION. 


The differential equation (3) is equivalent to the 
two equations 





dy _ Be oe 


dx a/ 1—2? 


Saas a (1—2?) + b (1—2?)9/? 
dx 


where g = sin ©. Assuming for yand z two power 
series * with undetermined coefficients, using the 
initial conditions x=0, y=0, zs = O and determining 
these coefficients so as to satisfy the above two 
equations, we deduce for y the expression 
2 4 
y= etd) i a Hl 5) 


“This method is suggested for approximate solutions of 
differential equations in RUNGE UND KONIG: Numeris- 


ches Rechnen. § 87. pp 321—23 
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Asa first approximation we can retain only the second 
term in (5), and the equation then illustrates how 
the assumption of non-uniform horizontal distribu- 
tion of g, affects the shape of the parabolic cable. 


I shall now indicate an alternative method of 
deriving equation (5) which is more illuminating. 
Equation (2), vez. 


d 
HG +45 (2) 


gives us the parabolic cable when we write @ in 
place of a V1 + p2 which is equivalent to putting 
p= 0 z¢e., treating the cable as straight. We can 
therefore obtain a better approximation by putting 








for p A ay, its valueas obtained for a parabolic 
cable. From § 1, the equation 


b x2 
y= : gives p = =bx and hence (2) can be 





approximated to 





anes 22 
Siege At +bx)t+ 5 

Te eds) b? x2 b* x’ : 

t.€.. — =ajylt — ere | 
Aint 2 3 as 2 


on expanding by the binomial theorem, since bx 
giving the slope at any point can be assumed, for all 
practical purposes, numerically less than unity, and 


the expansion being therefore valid. Retaining 
only two terms, this gives 
d*y ab? x2 





dx =(atb) + 


Integrating once and using dy — ifx =0 
dx 
dy 


—=(¢ + b)x + ab? x 
dx 6 
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Another integration, using y = 0 if x = 0 gives 


i (a + b) x? ab" x 5 
y 2 a" Toa: (5) 
which is the same * as equation (5). For use in 
practical calculations it is preferable to express this 
equation in terms of fand 7; but we have only one 
condition expressing that (5) passes through the 
points of support ( + 27/,, ¢)and_ this does not enable 
the determination of both aand 6 in terms of these 
quantities. Let us, however, write a=bn, so that n 
is the fraction of the cable weight to the horizontal 
load; equation (5) then reduces to 
y= b(n+1) x? nb" x 
2 24 


Putting x=1/2 & y=f, we have 


ris b (n+1)/ nb l 
8 384 
2 2 2 72 
l — ble bel (6) 


b 
Pom ie eal 
a er eee 48 





This equation gives a cubic in 6 from which we can 
find b in terms of n, f and 7 which might be supposed 
to be known quantities. We can solve this cubic 


equation approximately. 


2 
Writing =), eqn. (6) reduces to 


a 402 \ pe 

f-v= nd(i aie , a cubic ind 

For small values of n, a first approximation is 
8f ° 


fHd& b => 





* Gee W. H. Burr: ibid p..69 where a similar method is 
adopted in a different connection, 
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which is a well-known expression in a parabolic cable. 
To obtain a second approximation let us substitute 
this value ) — fin the right hand side and obtain 


) =s-nf (1+ $4)= a Or. -af 


. » os _. 4nfe 
4. Gey b= 72 ben ceR (7) 
8 8 
Ala; ee hal ley eee (3) 
dx 6 


6b being given from (7) in terms of n, f and J. 


Let us now calculate the length of are of the 
entire cable as modified by the above formule and 
compare the same with the original expression 


L=1 41+ mia 6% Yi 
Ey ok ie Py 


obtained in § 1. Denoting the modified length as 
L’ we have 


L/2 2 24 ; : 
Lig fT] 1+(# 3 oe being given from (8) 


= f+ ee pabie’ Her 


Expanding By the binomial theorem and retaining 
terms up to x‘, this reduces to 


va2 ff 4 ee { (at b)Aabt— 3(a+b) i, “Je 


an [+ (a+b)? | (a+b) § 4ab2—3 (a +b)’ F 
24 1920 
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on performing the integration indicated. Now, to 





substitute @= bn and b= p(i-« —2t in this ex- 


pression, we can write 
—64f7f, <4nf2\? 
pat ( nN 372 


=48(1 sedis <A OFT © 


<n ~“377 J» tetaining only first power of n, 











bee s_ 642 _ 16 nf? 
Similarly b°= 3 1—4n sah) 


Nope loz [ b? “s baeeatieee - PE +n) 


(3n°+3 n —n+)3] 


4 74 


1920 


Co 





=f! $14 2+ an) (3+2n)} 


° . 3 
neglecting terms in 12, 1’, etc. 


ee role a 14 aE (a- oa) | 
Bene 54" 1920 12 so) { 


Substituting in this the values of 6? and 6‘ written 
above, we get finally 


6 
) g f2'" 3aif4 | 148 nb pia fh 
tea ee epi Me La 8 








Hence L’ - L=1 = n f4/14 (neglecting the last term) 
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In practice it is found that all ratios of fit will 
be between .2 and .05 and the following table * 
gives values of L/l corresponding to different values 


of #/1. 





f/l Lil 
v2 1.09643 
yi75 1.07566 
15 1.05676 
A253 1.040104 
4 | 1.026027 
02> | 1.014797 
05 1.006627 





* Vide Burr : ibid. (p. 73) 
Siew 

l eke 7. 
gives us an indication of the error committed in 
treating go as a uniform horizontal load. Taking from 
the table in the Introduction the case of the Freiburg 
suspension bridge for which n=} (approximately) let 
us calculate the deviation for the case s/J=.2; we 
easily find 

- 7 F = 005262 

such that the corrected value would be 1.10169, 
showing that although the modification required is 
not great, it is of the same order as the difference 
between the values of L/z calculated as in the above 
table and the value calculated from the formula + 
exact for a parabolic cable, viz. 


1 
Lii= Ten (2! + Sin h2u ) 


The above expression 











+ Vide Burr: ibid. (p, 7) 
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§ 3. OTHER DIFFERENTIAL EQUATIONS 


In the previous investigation we have considered 
the weight of the cable alone and neglected that of 
the suspenders. Let g’ be the weight of the suspenders 
per square unit of vertical plane of the cable so that 
g’y is weight per linear horizontal foot of that part of 
the suspenders above the vertex of the cable. The 
part of the suspenders below the vertex forms a 
uniform load over the span and might beassumed to 
have been included in W,. The general differential 
equation becomes 








2 

T= aN 14 ph tb + ey 

; d ‘ 
ie, baa i452 +b +e (c=g'/H ) 


Making # and y the dependent and independent 
variables respectively, this can be reduced to the first 
order differential equation. 





pee V1 +p°+b+cy 


of which an exact solution is complicated. It can 

however be solved by approximate methods by 

substituting for y and 4» the expressions as given 
dx 


by the parabolic cable. 


2 f 
Putting y="sand px, equation (10) reduces to 


2 Pea Y, 2 
ao /; 1+b2 x2to+ Gus 
dx? 2 





=a ( 1+2E) ++ oes 


=(a+o)+ °F (b+¢) 
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Successive integrations give 


dy _ Peabo: bx. (ab+c) 
dx 6 

aa (at+b)x? , (abtc) x’ 

2 24 


which expresses the deviation from the parabolic 
form. 


I now consider, finally the case where the cable 
is not necessarily uniform. Let g, be the weight per 
ft. run at the vertex and g, at any point where 
inclination of the cable to the horizontal is 6, 
supposing that the normal section of the cable at 
any point is proportional to its total stress, we have 


& x=. SecO x SecO =g, Sec? ©. 
and the differential equation becomes 


dy _ 

——=a Se®O+b+cy (11) 
dx? 

Putting a2 = °F ¥and y= +f x? 


this can be written 


dzy_ 64 f2,4cf 
dx? plat b) + (zs re x2 


Se 2 
giving aes (a+b) +k x! 


F ee ol Ot fee acs 

where k= 12 (2 +e) 
An exact solution of equation (11) is a matter of 
some difficulty. 


I hope to return in a further paper to more 
detailed discussions of equations (10) and (11) and 
also discuss, if any, the bearing of the results of § 2 
on the problems of stiffening of the structure and the 
deflection of the bridge system, 
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